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i 1 ABSTRACT 

*-£h ■ We present CCD imaging observations of early-type galaxies with dark lanes obtained 

with the Southern African Large Telescope (SALT) during its performance- verification 
Q . phase. The observations were performed in six spectral bands that span the spectral 

S-^ ' range from the near-ultraviolet atmospheric cutoff to the near-infrared. We derive 

c/3 \ the extinction law by the extragalactic dust in the dark lanes in the spectral range 

l.ll/im -1 < A -1 < 2.94/«n -1 by fitting model galaxies to the unextinguished parts 
of the image, and subtracting from these the actual images. This procedure allows 
I , the derivation, with reasonably high signal-to-noise, of the extinction in each spectral 

J> ■ band we used for each resolution element of the image. We also introduce an alterna- 

■^j- ' tive method to derive the extinction values by comparing various colour-indices maps 

7— I , under the assumption of negligible intrinsic colour gradients in these galaxies. We than 

C^ - ■ compare the results obtained using these two methods. 

We compare the total-to-selective extinction derived for these galaxies with previ- 
. ously obtained results and with similar extinction values of Milky Way dust to derive 

f^) ' conclusions about the properties of extragalactic dust in different objects and condi- 

OO , tions. 

■ We find that the extinction curves run parallel to the Galactic extinction curve, 
^ ' which implies that the properties of dust in the extragalactic enviroment are similar 

• i— I , to those of the Milky Way, despite our original expectations. The ratio of the total V 

■ band extinction to the selective extinction between the V and B bands is derived for 

each galaxy with an average of 2.82 ±0.38, compared to a canonical value of 3.1 for the 
Milky Way. The similar values imply that galaxies with well-defined dark lanes have 
characteristic dust grain sizes similar to those of Galactic dust. We use total optical 
extinction values to estimate the dust mass for each galaxy, compare these with dust 
masses derived from IRAS measurements, and find them in the range 10 4 to 1O 7 M . 
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INTRODUCTION 

Since its existence was first inferred from the obscuration of 
starlight (Trumpler 1930), interstellar dust has been recog- 
nized as an important component of the interstellar medium 
(ISM) playing role in the evolution of galaxies, the forma- 
tion of stars and planetary systems, and possibly, the origins 
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of life. Studying the properties of the dust particles in in- 
terstellar space can help understand the nature of the dust 
grains themselves, the processes that govern their evolution 
and the way those grains affect the light of stars in their 
vicinity. The last property is mainly due to the interaction 
of dust grains with the electromagnetic radiation, i.e., due to 
the scattering and absorption of light, which depend strongly 
on the size distribution, structure and chemical composition 
of the grains. This attenuation of starlight is wavelength- 
dependent and is expressed as the wavelength dependence 
of the extinction, or the "extinction law". Since extinction 
tends to be greater in blue than in red, this is often referred 
to as "reddening". 
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The extinction law in the Milky Way (MW) has been 
studied extensively for many decades since it was first mea- 
sured for the optical region by Rudnick (1936) using the still- 
widely-used "pair match" . This method is based on the com- 
parison of the spectral energy distribution (SED) of pairs of 
stars of exactly the same spectral and luminosity classes. 
One of the stars is usually selected to be nearby thus is as- 
sumed to be clear of dust with a very small chance of being 
significantly extinguished. Any difference between the two 
SEDs is attributed to dust extinction and the global extinc- 
tion law is derived after comparing many pairs of stars. From 
the far-UV to the far-IR, the standard MW extinction law 
is generally adopted to be that derived by Savage & Mathis 
(1979; hereinafter S&M). However, there are significant dif- 
ferences between different lines of sight that cause large vari- 
ations. The S&M law was revised by Cardelli et al. (1988, 
1989) and Rv, the ratio of total to selective extinction, was 
proposed as the single-parameter value determining it and 
varying in different environments. While Rv is typically 3.1, 
values as low as 2.1 and as high as 5.6 have been measured 
(Valencic et al. 2004) . An explanation of the general extinc- 
tion law as a manifestation of the size distribution of dust 
grains was proposed by Greenberg & Chlewicky (1983, see 
their Fig. 1). Measuring the wavelength dependence of the 
extinction offers, therefore, a means of estimating the size 
of dust grains in the MW and in other galaxies. 

The extinction law in other galaxies has been measured 
using the pair method only for relatively nearby objects, 
such as the Large and Small Magellanic Clouds where indi- 
vidual stars can be resolved (e.g., Gordon et al. 2003), and 
significant differences from the MW extinction law were ob- 
served. As far as other galaxies are concerened, the above 
method is usually not practical due to the difficulty to re- 
solve individual stars, thus other methods must be used. 
Some of the methods suggested in the literature includes 
measuring the extinction by foreground galaxies of distant 
quasars (Ostman, Goobar & Mortsell 2008), determining 
the differential extinction of multiply-imaged quasars (Falco 
et al. 1999), determining the dust extinction of gamma-ray 
burst host galaxies (Li, Li & Wei 2007) and measuring the 
extinction law in starburst galaxies (Calzetti et al. 2000). 
However, in this article we adopt a different method, which 
studies the extinction law in early-type galaxies with dust 
lanes. 

The existence of early-type galaxies with dust lanes was 
recognized following the identification and study of five el- 
lipticals with minor-axis dust lanes by Bertola & Galleta 
(1978), and the publications of galaxy catalogs for objects 
with dust lanes by Hawarden et al. (1981), Ebneter & Balick 
(1985) and Veron-Cetty & Veron (1988). Moreover, many 
early-type galaxies show the presence of a multiphase ISM 
with interstellar dust as an important constituent. This is 
indicated by recent surveys revealing the presence of dust 
extinction in a large fraction of early-type galaxies (Van 
Dokkum & Franx 1995; Ferrari 1999; Tomita et al. 2000; 
Tran et al. 2001). The method is based on the assump- 
tion that the underlying galaxy, whether elliptical or lenticu- 
lar, has a smoothly-varying light distribution while the dust 
lanes are only local disturbances of the global brightness pat- 
tern. If the dust covering factor of the galaxy is relatively 
small, one can extract the dust-free model of the original 
underlying galaxy. Comparing the light distribution of the 



unextinguished galaxy with the actual observed one allows 
the derivation of extinction at each measured angular resolu- 
tion element of the target galaxy. This method was used for 
the first time to derive the properties of dust in the lenticu- 
lar galaxy N7070A (Brosch et al. 1985) . A discussion of the 
extragalactic extinction law in early-type galaxies, follow- 
ing a number of other studies, was summarized by Brosch 
(1988). Later attempts to derive extragalactic extinction law 
in early- type galaxies were by Brosch et al. (1990), Brosch 
and Loinger (1991) Goudfrooij et al. (1994), Sahu, Padney & 
Kembhavi (1998), Patil et al. (2002) and Patil et al. (2007). 

Brosch & Loinger (1991) used a simple elliptical- 
isophote fit to model the unextinguished part of the galaxy 
NGC 7625 and derived an extinction law rising faster at its 
blue end than the accepted MW standard. In latter publi- 
cations by Goudfrooij et al. (1994), Sahu et al. (1998), Patil 
et al. (2002) and Patil et al. (2007), the light distribution 
of the underlying galaxy was modeled after first masking off 
the dust lanes. This was done first for the reddest spectral 
band (usually the I band) by allowing a free-parameter fit 
to determine the optical center of the galaxy in a band that 
is least obscured by dust. This fit was repeated for the bluer 
bands keeping the galactic center fixed. This way, the ex- 
tinction at a specific location in the galaxy was determined 
for each wavelength used in the observations. The overall 
extinction laws were compared with the standard MW ex- 
tinction law. 

One overall result found in most studies of extragalactic 
dust is its similarity to the dust in the MW, at least in 
the properties revealed in the optical part of the spectrum. 
This is amazing, considering the wide variety of galaxies and, 
presumably, of physical conditions of the dust. Unless some 
cosmic plot is afoot, one would expect dust grains of different 
sizes and different chemical compositions in different objects 
- yet this does not appear to be the case. 

While a significant number of early-type galaxies have 
been investigated with this method (e.g., 26 objects in Patil 
et al. 2007 and 10 objects in Goudfrooij et al. 1994, none in 
common) , in most cases the measurements were restricted to 
bands longwards of B. As noted by Nandy (1984) from stud- 
ies of the Magellanic Clouds extinction, differences among 
extinction curves are better detected at shorter wavelength 
such as the U-band or the space ultraviolet. This paper aims 
to extend the optical study of the extinction law in a sam- 
ple of galaxies from the standard B, V, R & I bands to 
the near-UV by using two filters that split the U-band into 
two intermediate spectral regions. This allows the extrac- 
tion of more information about the extragalactic extinction 
law close to the blue atmospheric cutoff. The observations 
were performed with a very large telescope, which allowed 
a reasonable sensitivity even at the shortest spectral band 
used, despite the atmospheric extinction. 

This paper is organized as follows: §[l]gives a description 
of all the observations and data reduction. In §[2] we present 
our results, analyze them in § [3] and present our interpreta- 
tion in § [4] We discuss the results in § [5] and summarize our 
conclusions in § [6] 
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Table 1. Global parameters for galaxies in our sample. 



Object 


RA 
(J2000.0) 


DEC 
(J2000.0) 


Morph. 

(RC3) 


B° 


v Helio 

(km/s) 


Size 
(arcmin) 


NGC3283 


10:31:11 


-46:15:02 


SO 


12.50 


2889 


2.6x1.5 


NGC3497 


11:07:18 


-19:28:21 


SO 


13.03 


3707 


2.6x1.4 


AM1118-290 


11:20:53 


-29:24:09 


SO 


14.22 


8197 


1.1x0.7 


NGC4370 


12:24:55 


+07:26:38 


Sa 


13.48 


722 


1.4x0.7 


AM1352-333 


13:55:34 


-33:54:05 


SO 


14.30 


4577 


1.1x0.9 


NGC5626 


14:29:49 


-29:44:56 


SA0+ 


13.98 


6853 


1.2x0.9 


AM1459-722 


15:04:35 


-72:37:00 


S? 


14.05 


4493 


1.3x0.9 


NGC5799 


15:05:36 


-72:25:57 


SO 


13.32 


3099 


1.3x1.0 


NGC5898 


15:18:13 


-24:05:49 


E0 


12.49 


2209 


2.2x2.0 



Table 2. Observing log. 



Galaxy Exposure time (sec) point source FWHM (") 





Date 


B 


V 


R 


1 


i;2 


Ul 


B 


V 


R 


1 


U2 


Ul 


NGC3283 


16.03.2007 


20(7) 


20(7) 


10(7) 


10(7) 


120(7) 


120(7) 


1.82 


1.74 


1.84 


1.82 


2.09 


1.85 


NGC3497 


24.02.2007 


30(6) 


30(6) 


20(5) 


20(5) 


120(5) 


120(4) 


2.38 


2.30 


2.41 


2.18 


2.41 


2.38 


AM1118-290 


26.02.2007 


30(5) 


30(5) 


20(5) 


20(5) 


120(4) 


120(4) 


2.02 


1.74 


1.59 


1.67 


1.92 


2.11 


NGC4370 


23.02.2007 


30(8) 


30(6) 


20(8) 


20(8) 


120(8) 


120(14) 


2.27 


2.37 


2.04 


2.13 


2.35 


2.35 


AM1352-333 


26.02.2007 


30(5) 


30(5) 


20(5) 


20(5) 


120(4) 


120(9) 


1.96 


1.88 


1.88 


1.85 


1.79 


2.07 


NGC5626 


17.05.2007 


30(9) 


30(10) 


20(8) 


20(7) 


120(6) 


120(6) 


2.49 


2.41 


2.27 


2.60 


2.55 


2.66 


AM1459-722 


16.03.2007 


30(9) 


30(10) 


20(9) 


20(10) 


120(10) 


120(9) 


1.90 


1.96 


1.88 


1.90 


1.51 


1.57 


NGC5799 


23.02.2007 


20(8) 


20(8) 


10(8) 


10(8) 


120(8) 


120(8) 


2.32 


2.27 


2.21 


2.18 


2.24 


2.60 


NGC5898 


21.02.2007 


6(6) 


18(6) 


25(6) 


45(6) 


76(6) 


92(6) 


1.88 


1.9 


1.65 


1.73 


1.73 


1.89 



Note: The number in parenthesis represents the number of individual frames combined to produce the final galaxy image in 
each band. The point source FWHM represents the typical FWHM of stars near each target combining seeing and image quality 
effects. 



1 OBSERVATIONS AND DATA REDUCTION 

The Southern African Large Telescope (SALT) was de- 
scribed by Buckley, Swart & Meiring (2006), and its CCD 
camera (SALTICAM) was described by O'Donoghue et al. 
(2006). We used the SALT and SALTICAM to observe dark- 
lane early- type galaxies. 

It may seem superfluous to observe bright galaxies with 
an (effectively) nine-meter telescope since smaller telescopes 
served admirably this purpose in the past. These observa- 
tions, however, were obtained during the Performance Ver- 
ification (PV) phase of the SALT telescope when only the 
SALTICAM was available, and even this with the instru- 
ment before its upgrade to auto-guiding. The telescope was 
still affected by the limited image quality, which restricted 
the size of the workable field to some three arcmin instead of 
the nominal eight arcmin region. The specific advantage in 
performing this project with SALT and SALTICAM is the 
high efficiency of this combination in the near-ultraviolet, 
and the availability of the Ul and U2 filters (see below) 
among the SALTICAM filter complement. 

The SALTICAM is a CCD mosaic of two 2048x4102 
pixels CCDs and therefore, for consistency, program objects 
were imaged with the same CCD throughout the observa- 
tions. The full SALTICAM image is about 9'.6 by 9'.6 with 
pixels of ~0".28 (after binning on-chip by a factor of two), 
but the nominal science field is eight arcmin. 

Observations of nine dust lane early-type galaxies were 
performed in service mode from February to May 2007 with 
SALT, whereas only for eight objects we derived the extinc- 




Wavelengthfnm] 

Figure 1. SALTICAM short wavelength filters: 340/45 nm (Ul- 
dashed line) and 380/40 nm (U2- solid line) transmission. 



tion values (see § [2J . The nine galaxies are listed in Table [T] 
with their coordinates, morphological classification and op- 
tical size as taken from the RC3 catalog (de Vaucouleurs et 
al. 1992) or from LEDA | |http://leda.univ4yonllr) . CCD 
imaging observations were performed with SALTICAM us- 
ing the standard B, V, R & I filters, as well as two short- 
wavelength interference filters, Ul and U2, with transmis- 
sion bands peaking at 340 nm (FWHM 35 nm) and 380 nm 
(FWHM 40 nm) respectively, as shown in Fig.Q] The plotted 
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transmission profile is logarithmic, to emphasize the reason- 
ably good blocking of out-of-band wavelengths. Ul has its 
strongest out-of-band transmission band above 1 /im where 
the SALTICAM CCD already has a fairly low quantum ef- 
ficiency. U2 is slightly more affected by a red "leak" given 
its prominent transmissions at ~500 nm and at ~700 nm. 
However, the first out-of-band peak at 480 nm is at less than 
1% level, and the second one, a blue leak near 310nm, is at 
~ 0.1% level. 

The objects were selected from catalogs for objects with 
dust lanes by Hawarden et al. (1981), Ebneter & Balick 
(1985) and Veron-Cetty & Veron (1988). The observed ob- 
jects were chosen according to visibility conditions and avail- 
ability of observing nights requiring that the galaxies would 
fit in the good image quality field of SALTICAM (~ three 
arcmin), would be in the sky region accessible to SALT 
at the date of observation, and would not present reduc- 
tion complications such as dust lanes with bright edges that 
might signify internally-produced light within the dust lanes. 
Galaxies from Patil et al. (2007) and from Goudfrooij et al. 
(1994) need the full observation set even though they were 
already observed by these authors in some of the bands. This 
is because the images to be used for the derivation of the 
extinction should all have the same seeing size. This is en- 
sured by obtaining the six images quasi-simultaneously, with 
the same camera and telescope. The observations consist of 
short-duration integrations while the filter wheel cycles be- 
tween the six filters, thus ensuring quasi-similiar seeing size 
for all bands observed in a single filter cycle. Exposure times 
are determined to be sufficiently short to prevent image trail- 
ing when operating in the SALTICAM unguided mode (but 
with SALT tracking) or overexposure of the brightest part 
of each galaxy. The six-image cycles were repeated while the 
pointing was slightly shifted to allow for a dithering effect. 
The cycles were then accumulated to reach better exposure 
depth. Exposure times in seconds for B, V, R, I, U2 & Ul 
bands, as well as the typical F WHM of stars near target in 
the different bands, are listed in Table [2] The numbers in 
parenthesis represent the number of exposures acquired in 
a given band. Since each image consists of a combination 
of a number of short-duration exposures, the FWHM value 
represents the value to which the different images in each 
band were convolved. The size of the resolution elements is 
a result of both seeing and image quality (IQ) effects. 

The SALT pipeline includes bias subtraction, overscan 
subtraction and "cross-talk" correction, while other stan- 
dard preprocessing steps were done using standard tasks 
within Such steps include geometric alignment of 

multiple frames taken in each filter by measuring centeroids 
of several common stars in the galaxy frames. The frames 
were then combined with median scaling to improve the 
S/N ratio. This alignment procedure involves IRAF tasks 
for scaling, translation and rotation of the images, so that 
a small amount of smearing is introduced affecting the ac- 
curacy to be better than half a pixel. Median combination 
was also useful in removing cosmic rays events while the 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories (NOAO), which is operated by the Association of Uni- 
versities, Inc. (AURA) under co-operative agreement with the 
National Science Foundation 



CCD hot pixels were removed with the CCDMASK task in 
IRAF using an appropriate mask. We emphasize that no 
Hatfield correction was made during this reduction process 
since at the time of observations the telescope did not yet 
have a moving baffle to simulate the effect of its continuously 
changing pupil on the flatfield during observations. 

Since no flatfield correction was applied, the sky back- 
ground cannot be estimated naivly by measuring at various 
locations in the frame away from the galaxy, or by fitting a 
tilted plane on the basis of small areas far from the galaxy 
images. In order to determine the sky background we took 
advantage of its statistical nature; we first created a count- 
frequency histogram for the galaxy surroundings (that is, in 
a 3' radius from the galaxy center) and then fitted a Gaus- 
sian distribution curve to the peak of the histogram, where 
the most frequent count value was assumed to represent the 
sky background. This value was subtracted from the galaxy 
frame. Due to the nature of the unflatfielded frames we 
avoided combining frames from different observing nights. 
Furthermore, only observations with similar seeing values 
within the same observing night were selected to produce 
the combined image after convolution with a Gaussian to 
match the image seeing. 



2 RESULTS 

As mentioned above, the observations reported here were 
obtained during the PV phase of SALT in a "shared-risk" 
mode. This implied that the telescope was not optimized 
for imaging, the IQ problem that troubled the telescope at 
that stage of commissioning was not yet solved and the mir- 
ror segments were not freshly aluminized. The implications 
are that the images outside a ~ 3' circle had reduced res- 
olution, and sometimes the seeing size degraded during the 
observation yielding some unusable images. 

We present contour maps corresponding to the various 
filters in use for each galaxy. Figs.[6]to[l4]show filled contour 
maps for each of the sample galaxies. The contour maps are 
created from the reduced and combined final images and are 
each plotted with 0.5 magnitude steps. As evident from the 
contour maps, the Ul and U2 images are significally fainter 
than the optical band images, producing at the dark lanes a 
S/N ratio of ~ 10 and ~ 25 per resolution element, respec- 
tively. This, while the characteristic S/N ratio per resolution 
element for the B, V, R & I images may reach values of ~ 100 
and even higher. Since the images are not calibrated, the es- 
timation of the contour level where the surface magnitude 
reaches to 25 mag/arcsec 2 is based on the semi-major axes 
measurements from the RC3 catalog (see de Vaucouleurs et 
al. 1992). The corresponding B-image contours are labelled, 
allowing for a comparison between the S/N ratio in different 
filters for each galaxy. Galaxies where the major-axis is too 
large to fit the 112"xll2" contour maps are not labelled. 
Some of the contour maps, mainly for the near-UV bands, 
show low S/N chevron vintage patterns due to both electron- 
ics and the unflatfielded nature of the images. In some cases, 
stellar images in Figs. [6] to [14] appear somewhat distorted 
mainly towards the edges, probably caused by the IQ prob- 
lems mentioned earlier. We also present B-I colour-index 
maps of the sample galaxies in Fig. 1151 These images show 
the morphological differences between the suspected dust 
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lanes, presenting concentric rings of dust (e.g., NGC3497 
and NGC5626), warped dust disks (e.g., NGC4370) and an 
orthogonal set of dust lanes (AM1459-722). 



3 ANALYSIS 

In order to estimate the dust extinction in galaxies one can 
compare the light distribution in the observed galaxy with 
its dust free model. Since the target galaxies are early-type, 
one could assume that they are spheroids and use a de Vau- 
couleurs (1948) light distribution 



H (r) = ^ + 8.325 



1/4 



- 1 



(1) 



or a more flexible fit using a Sersic light distribution for a 
bulge which would fit an elliptical galaxy 



Hb (x,y) = fie + k n 



l/n 



(2) 



where r&= 



-.[\x 



c+2 



+ 



, qb is the axial ratio of the 



1 lb 1 J 

bulge and a pure ellipse has c=0. A more realistic fit assumes 
that the galaxy is a superposition of a bulge following the 
Sersic law and an exponential disk, whose light distribution 
is given by 



/id (x, y) = fi + 1-086 ( T -j- 



(3) 



where r d = [j2;| c+2 + | ^ | c+2 j <c+2) , q d is the axial ratio of the 
disk and a pure ellipse has c=0. Such an assumption would 
fit a lenticular galaxy. 

One method of fitting a smooth galaxy model to the im- 
ages is to use the "standard" model mentioned above for the 
entire galaxy image, which may be preceded by the masking 
of parts of the image where obvious dust is present. It has 
also been shown that constructing a dust-free model for an 
early-type galaxy, which has a fairly smooth and symmetric 
light distribution with respect to its nucleus, can be easily 
done by fitting ellipses to the isophotes of the observed im- 
age. Although this method has been used succesfully in a 
number of cases (Brosch and Loinger 1991; Goudfrooij et 
al. 1994; Sahu et al. 1998 and Patil et al. 2007), it seems to 
be less effective in our case. First, in most cases, the mask- 
ing of the entire dusty region in the elliptical fits leads to a 
percentage of masked pixels that is too large to be handled 
by the ellipse- fitting routine in IRAF. Second, it seems that 
fitting ellipses for the images in the near-UV produces poor 
results mainly due to the characteristic relatively low S/N 
ratio per resolution element (which can even be as low as 5) 
following the background subtraction. For these reasons we 
selected a different fitting method that is described below. 

3.1 Extinction maps 

The appearance of the dust lane, regardless of its orientation 
with respect to the stellar body, can be described in several 
ways (Bertola 1987). The dust lane may lie exactly along 
one of the axes of the stellar body, presumably representing 
a disk seen edge-on, or show a slight curvature, suggesting a 
disk seen almost edge-on. Dust lanes may also appear as full 



rings or may appear to be warped at their outer parts where 
the minor-axis dust lane bends toward the major axis in the 
outer regions. In a few cases, a set of multiple parallel dust 
lanes is present, suggesting a system of coplanar rings seen 
at an angle. Considering the dust morphology, one should 
also consider the extent of the dust lane with respect to the 
galaxy, i.e., whether it is confined to the innermost regions of 
the galaxy or whether it follows the luminous stellar material 
down to its detection limit. 

While the dust lanes do not always appear to be sym- 
metric under reflection through the nucleus, we may assume 
that the light distribution of the underlying galaxy does hold 
such symmetry properties in order to extract the extinction 
maps. The extinction at each location within a resolution 
element may be measured using 



Ax 



-2.5log 



I\,ob 



(4) 



which under reflection symmetry may be expressed as 

A x = -2.5^ [loglx (xi,y s ) - logh (-*», ~Vj)] (5) 

ij 

where Ax represents the amount of extinction in a particu- 
lar band on a magnitude scale, and 7a (x, y) represents the 
average intensity of a rectangular seeing-size box at location 
(x, y) away from the nucleus. Such rectangular boxes were 
translated over the dust-occupied regions of each galaxy 
with no overlap, though the nuclear regions (radius^5") 
were excluded to avoid seeing-related effects. We also as- 
sume the Galactic extinction in the line of sight is uniform 
over each galaxy, thus no correction for Galactic extinction 
is necessary. 

However, the reflection symmetry does not necessar- 
ily hold when major-axis dust-lane galaxies are concerned. 
Three of our sample galaxies (NGC3283, NGC3497 and 
NGC4370) were studied by Bertola et al. (1988) as part of 
their study of major-axis dust-lane ellipticals. While from 
the B-I colour-index maps it appears that the dust lane in 
NGC3497 does not lie exactly along the major axis (see also 
Patil et al. 2007), we do identify a major-axis dust lane in 
NGC4370 and NGC3283. We should note here that there is 
a tendency in the usual morphology classification schemes 
to call SO any galaxy where a stellar body is crossed by a 
dust lane, regardless of its orientation. 

In order to extract an extinction map for each of these 
two galaxies we used a two-dimensional bulge/disk decom- 
position analysis code (BUDDA; Souza, Gadotti & dos An- 
jos 2004) to derive the galaxies' global parameters and to 
build 2D models for the galaxies. In order to determine ini- 
tial values for the galactic center, position angle, ellipticity 
and ellipse index (i.e., the galaxy boxiness), we fitted ellipti- 
cal isophotes to the outer regions of the galaxies, which are 
least affected by dust. Initial values for the effective radius 
and effective intensity (see eq. [TJ were obtained by fitting 
a de Vaucouleurs profile to the minor axis, which was first 
rotated to lie parallel to one of the x-y axes for convenience. 
Excluding the inner regions occupied by dust, this profile 
proved to fit well the light distribution along the minor axis, 
supporting the assumption that these galaxies are indeed el- 
lipticals. This fit was initialy performed for the optical band 
images in an iterative mode, meaning that dust regions were 
excluded from each fit according to the residual image based 
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on the previous fit, until revealing what seems to be the full 
extent of the dust lanes. For consistency, we compared the 
global parameters obtained for each band, including the lo- 
cation of the galactic center, the position angle, the ellip- 
ticity, the ellipse index and also the effective radius. These 
values were found to match and were later used to fit models 
for the Ul and U2 images, where only the effective intensity 
was left as a free parameter due to the relatively low S/N 
ratio of these images. 




3.2 Extinction curves 



As explained in § 13.11 total extinction values were measured 
for different dust-extinguished regions in each galaxy. Those 
values were later used to derive the extinction values for each 
galaxy by fitting a linear regression between the total extinc- 
tion A\ and the selective extinction E(B-~V)=Ab — Ay using 
weighted least-squares with a fixed zero extinction point. 
The best-fitting slopes were used as the extinction values 

R x (= E (gtv) ) f° r eacn g a l ax y an d were compared with the 
Galactic values. We also used this method to fit a linear re- 
gression between the total extinction at wavelength A (Ax) 
and the optical extinction Ay. Extinction values normal- 
ized to the optical extinction were obtained by averaging the 
best-fitted slope of A\ versus Ay and the reciprocal slope of 
Ay versus A\. Extinction values obtained using the latter 
normalization seem to produce less noisy profiles, probably 
as a result of the larger relative errors of E(B-V) compared 
to those of Ay, where the B image of an early- type galaxy 
is usually significally fainter than the V image. 

The derivation of the dust extinction law using the pro- 
cedure described above may be affected by the location of 
the dust within the galaxy. While in our Galaxy we assume 
all the dust is between us and the light source, this may not 
be the case for dust lanes or dust patches in galaxies where 
dust is located within the stellar body. This phenomenon can 
be illustrated using two simple models developed by Walter- 
bos & Kennicutt (1988) for two idealized cases. In one case, 
an optically and geometrically thin dust lane is embedded in 
the galaxy, thus the observed extinction A\ can be expressed 
as: 



Ay 



-2.5 log [p+(l -p)e 



(6) 



where p is the fraction of light emitted at wavelength A that 
originates between us and the dust, and rx is the optical 
depth of the dust expressed in units of tv, the optical depth 
in the V-band, using the Galactic reddening law to make 
the transformation. In the other simplified case, the dust 
and stars are homogeneously mixed in the dust lane, thus 
A\ can be expressed as: 



A x = -2.5 log 



1 - e" 
rx 



(7) 



where tx is computed on the same basis as in eq. [6] Con- 
sidering eq. [S] and eq. [7J these two effects seem to increase 
the slope of the observed extinction curve and the observed 
extinction values with respect to their intrinsic values in the 
optical part of the spectrum (see, for example, Brosch et al. 
1990). This deviation from the intrinsic extinction curve is 
more significant for larger optical depths. 




Figure 2. Effects of foreground light on the observed extinction. 
The dashed line is the typical S&M relation while the solid lines 
are for p = 0.05, p = 0.2 and p = 0.5 values. Optical depth values 
run between ry = and ry = 2. 



It is also interesting to study the influence of the dust lo- 
cation on the observed relation between Ax and E(B-V) for 
the case of varying optical depths, as may be found in differ- 
ent regions occupied by dust. The relation Ax ~ 1.086 tx is 
valid as long as we assume the dust is a foreground screen. 
Since for a given wavelength and a given dust grain dis- 
tribution the optical depth may only change with varying 
column density and column length (see also § Q , this rela- 
tion implies an expected linear correlation between Ax and 
E(B-V) where the slope represents the extinction value Rx- 
However, even relatively small amounts of light originating 
between us and the the dust lane may significantly affect the 
observed Ax to E(B-V) relation, as illustrated in Fig. [5] and 
Fig. [3] Therefore, one should be careful when fitting a linear 
regression to obtain an extinction value, where in the case 
presented above this linear relation may not be valid even as 
a first approximation. In addition to the Ax versus E(B-V) 
relation, which appears non-linear for relatively small val- 
ues of p and small optical depths, Fig.[5]also shows a similar 
effect on the Ax versus Ay relation, although as the figure 
shows, the deviation from linearity becomes significant only 
for high p values and very high optical depths. 

As is apparent from Fig. 1151 the spatial distribution of 
the dust is clearly not that of a foreground screen between 
us and the light sources. We may use a representative value 
of p = 0.15, estimated by Knapen et al. (1991) for the Som- 
brero galaxy, as a simple test case to estimate the effect 
of a dust lane located between us and the light source on 
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Figure 3. Effect of homogeneously mixed dust and stars on the 
observed extinction. The dashed line is the typical S&M relation 
while the solid line is for the embedded dust. Optical depth values 
run between ry = and ry = 2. 

the derived extinction law. Fitting a synthetic curve with 
p = 0.15, illustrated in Fig. [5] with a linear trend shows a 
very significant deviation from the S&M relation. However, 
an attempt to fit the sample galaxies under the assumption 
of p — 0.15 in order to obtain the optical depth in each dusty 
location does not show a significant deviation from the re- 
sults obtained with p — 0, with the R\ values only slightly 
decreasing. This may be explained if the optical depth in 
the dusty locations is relatively thin, i.e., ry < 1 as typical 
for our sample galaxies, where the deviation from the "real" 
linear relation may not be significant and within the error 
bars. We also note that due to the spatial resolution, only 
about a dozen dusty locations were measured within each 
galaxy, while it is possible that the dust geometry and pat- 
tern may change throughout the galaxy. Therefore, better 
spatial resolution is required in order to trace the possible 
variations of the dust lane location within the galaxy. One 
may reproduce the "real" linear relation from the observed 
one by fitting for different p values. However, considering 
a more complex dust geometry, this task would require the 
use of more complex radiative transfer models. 

We should also consider the possibility that part of the 
light from the dust lanes is due to the forward-scattering 
of light from stars in the galaxy. It is well-known that the 
amount of scattered light increases strongly towards shorter 
wavelength where the dust albedo is higher. Therefore, the 
scattering of light artificially increases the observed and 
subsequently "flattens" the extinction curve. Bruzual, Ma- 



gris & Calvet (1988) have treated the complex case of mul- 
tiple scattering in disk galaxies and demonstrated that, for 
an almost edge-on galaxy, Ry can reach values as high as 5 
for a reasonable amount of dust (also see Goudfrooij & de 
Jong 1995). Emsellem (1995) suggested that the scattering 
of light by dust particles plays a crucial role on the observed 
extinctions of the highly inclined Sombrero galaxy; this is by 
acting as an additional source of radiation and thus reduc- 
ing the observed attenuation by a factor of 2. He concluded 
that the relatively high extinction values of the Sombrero 
galaxy are due to the saturation of the observed attenua- 
tion for Av > 3 mag. Using Monte Carlo techniques for 
radiative-transfer calculations, Witt, Thronson & Capuano 
(1992) and Baes, Dejonghe & Davies (2005) showed that 
the effects of scattering can be important, even for small 
optical depths, and that neglecting or approximating the ef- 
fects may lead to serious extinction errors. These authors 
also showed that the forward scattering strongly depends 
on the geometry of the dust and star clouds. For a uniform 
sphere of dust within a larger uniform stellar distribution, 
Witt et al. (1992) showed that only a few percent of the 
unextinguished light is forward-scattered for a total extinc- 
tion typical for our sample galaxies (i.e., ~0.5%). Since the 
measured extinction values are not significally larger than 
the MW extinction values, and the extinction curves do not 
show any flattening towards the near-UV, it is fair to as- 
sume that the effect of scattering is negligible compared to 
the effect of the grain size, at least for those galaxies which 
are not highly inclined. 

We derived extinction curves and extinction values for 
eight of our nine sample galaxies, with the exception of 
NGC5898. The obtained R\ extinction values are listed in 
Table [3] along with the canonical extinction values for the 
B, V, R & I bands taken from S&M for comparison. The 
extinction curves are plotted in Fig. U and [5] For the case of 
the faint arc-like dust lane of NGC5898 (see B-I colour index 
map in Fig. 11 5 P , the measured V-band extinction per reso- 
lution element is rather small with values up to 0.1 mag. In 
addition, these galaxy images are relatively faint due to short 
exposure times (see Table [2} and therefore the extinction is 
poorly determined. For the case of NGC3283 and AMI 118- 
290, the B-I colour-index maps show prominent dust lanes 
along the galaxies major axes. These galaxies show noisy 
extinction plots in all bands, thus extinction values are dif- 
ficult to determine using only a linear regression fit, as has 
been suggested before. 

For the case of major-axis dust lanes it is possible that 
the dust lanes are actually dust disks seen edge on, therefore 
we should not expect the relation A\ = 1.086 t\ to remain 
valid. A 2D decomposition analysis of AM1118-290 shows a 
stellar disk that may be responsible for significant forward 
scattering from the major-axis dust lane. While NGC3283 
shows a simple de Vaucouleurs profile with a boxy struc- 
ture, as also indicated by Bertola et al. (1988), the same au- 
thors observed radio jets from the galaxy center, implying 
the presence of an active central source that could produce a 
significant forward-scattering from the major-axis dust lane. 

We should note that the derived Ry values for 
NGC3283 and AM1118-290 are ~2.42 and ~4.05, respec- 
tively. The first may imply that forward scattering, which 
tends to increase the extinction values, may not be dominant 
in the case of NGC3283. The extinction curves for these two 
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Table 3. R\ values 



Object 


Rui 


Rl72 


Rb 


Rv 


Rr 


Ri 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC3283 


4.14 ±0.03 


3.85 ±0.02 


3.42 ±0.01 


2.42 ±0.01 


1.98 ±0.01 


1.96 ± 0.01 


NGC3497 


4.41 ± 0.29 


4.10 ±0.19 


3.68 ±0.16 


2.68 ± 0.11 


2.11 ±0.09 


1.23 ± 0.05 


NGC4370 


5.07 ±0.09 


4.44 ±0.03 


3.89 ±0.02 


2.86 ± 0.01 


1.98 ±0.01 


1.05 ± 0.01 


AM1352-333 


4.57 ±0.32 


4.15 ±0.13 


3.61 ±0.07 


2.62 ± 0.02 


1.72 ±0.02 


0.98 ± 0.02 


AM1118-290 


7.17 ±0.46 


5.67 ±0.11 


5.02 ±0.06 


4.05 ± 0.03 


3.12 ±0.03 


2.26 ± 0.03 


NGC5626 


4.54 ±0.23 


3.93 ±0.09 


3.88 ±0.05 


2.92 ± 0.02 


2.23 ±0.02 


1.46 ± 0.02 


AM1459-722 


4.77 ±0.17 


4.46 ±0.05 


3.91 ±0.03 


3.08 ± 0.02 


2.32 ±0.02 


1.71 ± 0.01 


NGC5799 


3.82 ± 0.11 


3.43 ± 0.04 


2.89 ±0.03 


1.90 ±0.01 


1.24 ±0.02 


0.67 ± 0.04 


MW galaxy 


4.9 


4.56 


4.10 


3.10 


2.32 


1.5 



galaxies are plotted in Fig. [5] We also note that for the case 
of NGC4370, another major-axis dust lane elliptical, there 

were no signs of a luminous central body in our 2D decom- 

A\ 

position analysis and although the — — values measured in 

Av 

various dusty regions produce somewhat noisy plots, we ob- 
tain extinction values very similar to the Galactic extinction 
values. 



colour values throughout the galaxies and redden where the 
dark lanes appear. As before, we fit a linear regression be- 
tween the A-V and the B-V values, where the slope repre- 
sents the extinction value gig and the intercept with 

the A-V axis depends on the constant colour of the dust- 
free galaxy in the different bands. This means that if one 
can measure these colour values in different dusty regions, 
then one does not need to specifically measure the constant 
colours, but only to assume that such values do exist. 



3.3 An alternative method - comparing 
colour-index maps 

We suggest here an alternative method for deriving extinc- 
tion values, based on the comparison of colour-index maps 
of dusty early-type galaxies. As the traditional method de- 
scribed earlier, this alternative method also relies on the 
assumption of smoothly varying nature of the light distribu- 
tion in early-type galaxies, in which the dust lanes or dust 
patches are only local disturbances. Considering the de Vau- 
couleurs and Sersic profiles (eq. [1] and eq. [2} for early-type 
galaxies, one may assume that the brightness variation in 
different locations in the galaxy with respect to its central 
brightness is wavelength-independent, i.e., is similar for all 
bands. As a result, all colour-index maps of such a smooth 
early-type galaxy should exhibit constant values for the en- 
tire galaxy. The hidden assumption here is that the stellar 
population mix does not vary much across the galaxy; only 
the column density of stars presented to the line of sight 
changes with galacto-centric distance. This is justified by 
the shallow colour gradients measured for normal early-type 
galaxies (e.g., La Barbera et al. 2008) where the small colour 
gradients are explained as a metallicity effect. 

If dust is present, and if this dust produces wavelength- 
dependent extinction in the optical, a colour-index map, e.g., 
a B-V map, is expected to redden at the dusty locations 
with respect to the constant colour of the dust free regions. 
This reddening may be linked to the difference in extinctions 
Ab — Ay - Therefore, by comparing different colour values of 
dusty regions one may obtain certain extinction values and 
study the dust properties. 

We compared the various A-V maps with B-V maps for 
each galaxy by measuring the colour values in different re- 
gions in the galaxies. As expected, the maps show constant 



This method of deriving extinction values has the dis- 
advantage of being somewhat less informative in comparison 
with the traditional method outlined in § 13.21 Using this al- 
ternative method one cannot measure the extinction in each 
band nor to derive R\ values. The latter may be of impor- 
tance if one wishes to examine whether the target galaxy's 
extinction curve lies above or below the MW extinction 
curve (see also §[4]). Nevertheless, comparing the y\ ex ~ 

tinction values of different galaxies with the Galactic extinc- 
tion values we may examine whether these extinction curves 
run parallel to the Galactic extinction curve or whether they 
deviate from this linear trend, implying different extinction 
laws. Despite its pitfalls, the alternative method has the ad- 
vantage of being very easy-to-use and straightforward, with- 
out requiring any model fitting but only assuming that such 
a smooth surface brightness distribution and colour model 
do exist. Moreover, while the traditional method relies on 
the dust-free regions of the galaxy in order to model the en- 
tire underlying galaxy, the alternative method may be used 
without such requirements, relying solely on the dusty re- 
gions themselves. This means that the alternative method 
may be used even in cases where the dust obscures large frac- 
tions of the underlying galaxy, or in cases where the model- 
ing of the underlying galaxy is not simple. In addition, the 
alternative method may be favoured when the image reso- 
lution is poor and there are not sufficient datapoints to be 
fitted by a linear regression. In this case, one may have to de- 
rive the colour of the galaxy by measuring dust-free regions. 
We plotted the extinction values obtained using the alterna- 
tive method along with the extinction values obtained using 
the traditional method in Fig. [4] As these plots show, the 
alternative method produces extinction values with larger 
errors compared with those produced using the traditional 
method. 
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Figure 4. Extinction curves for the program galaxies (solid lines) along with the canonical curve for the galaxy 
(dashed lines) for comparison. The upper curves for each galaxy represent R\ values whereas the lower ones represent 
A\? 

values. Extinction values derived from previous papers are from Goudfrooij et al. (1994) and Patil et al. 

Ab — Av 

(2007), and are represented by crosses and open circles respectively. The error bars are la errors. 



4 INTERPRETATION OF DUST PROPERTIES 



4.1 The derivation of an extinction law 



We adopt the Mathis, Rumpl & Nordsieck (1977) dust grain 
model succesfully used in reproducing the MW extinction 
curve in the optical-UV part of the spectrum. By varying the 
grain size parameters and the abundance ratio of graphite 
and silicates we reproduce the galactic extinction curve and 
determine the dust grain characteristic size. We find the av- 
erage characteristic size in our sample galaxies to be similar 
to the characteristic MW dust grain size. We also estimate 
the dust mass in the sample galaxies from optical extinction 
measurements and find them to be significantly smaller than 
the dust mass derived using far-IR data. 



Assuming that the chemical composition of the extragalac- 
tic dust grains is uniform throughout each galaxy and is 
similar to that of the dust in our Galaxy, one can compute 
the relative grain size and the total dust mass in the sample 
galaxies by using available models for the chemical compo- 
sition and the shape of the dust grains. The term "interstel- 
lar dust" refers to different materials with rather different 
properties and there are significant differences in extinction 
between low-density diffuse dust and dust in the inner and 
outer parts of molecular clouds. Therefore, for the case of 
dust lane galaxies, we can only refer to a "mean" or typical 
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Figure 5. Extinction curves for NGC3283 and AM1118-290 (solid lines) along with the canonical curve for the 
galaxy (dashed lines) for comparison. These galaxies show noisy extinction plots, therefore the derivation of the 
extinction values is suspect. The upper curves for each galaxy represent R\ values whereas the lower ones represent 

value, The error bars are terrors. 
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extinction law, assumed to represent the large-scale diffuse 
dust. 

Since the total extinction linearly depends on the ex- 
tinction cross-section, we may express R\ as: 

Rx = C °*<W , ( 8 ) 

where C ex t represents the extinction cross-section of grains 
in a volume unit. The mean cross-section for spherical par- 
ticles with a size distribution of rii (a) da, where i indicates 
a certain component and a represents a certain grain size, 
can be written as: 

Si f a + Qext (a, A) ira 2 rii (a) da 
(C ext (A)) = (Q ext (a, A) vra ) = " Ei/;+Ma)da 

(9) 

where a_ and a+ represent the lower and upper cutoffs 
of the size distribution, respectively; A is the wavelength; 
Qext (a, A) is the ratio of the extinction cross-section to the 
geometrical cross-section and is denned as the extinction ef- 
ficiency. 

Several grain models, invoking different size distribu- 
tion functions and mixtures of chemical compositions, are 
proposed in the literature. Since we restrict ourselves to the 
extinction in the optical part of the spectrum, it is possible 
to use a simple model that can account for the MW optical 
extinction curve. We adopt the Mathis, Rumpl & Nordsieck 
(1977; hereinafter MRN) two-component model consisting of 
individual spherical silicate and graphite grains with an ad- 
equate mixture of sizes. The MRN model assumes uncoated 
refractory particles having a power-law size distribution of 

n(a) da = uhA{ a -3,5 da, (10) 

where a represents the grain size, uh is the hydrogen number 
density and Ai is the abundance of component i. The size 
distribution is truncated at the upper end at a+ = 0.25/im, 
with the lower end of sizes extending downward to a_ = 
0.005/im. The MRN model was also adopted by Draine & 
Lee (1984; hereinafter DL84) who succesfully used the Mie 



scattering theory for homogeneous spherical particles in or- 
der to reproduce the MW extinction curve in the optical- 
UV part of the spectrum. DL84 adopt a graphite-to-silicate 

abundance ratio of —r- ~ 1.12 to optimize the fit of their 
A s 

synthetic extinction curve to the observational data. We em- 
phasize that this model consists only of spherical particles, 
while observations of interstellar polariztion indicate that 
an appreciable fraction of the grains must be nonspherical. 
Since we do not deal here with polarization, we disregard 
non-spherical grains. 

Using eq. JS} and eq. (O, one can calculate the 
theoretical extinction values, where the extinction effi- 
ciency for specific wavelength, grain size and refractive 
index is calculated using the standard Mie theory (Mie 
1908; Debye 1909) given the dielectric functions for 
graphite and 'astronomical silicate' derived by DL84 
and published in tabular form (Draine 1985; see also 

: / /www.astro. princeton.edu/~draine /dust /dust .diel.html ) . 
The method for calculating the extinction values is also 
outlined by Steenman & The (1989). 

4.2 Dust characteristic size 

We now estimate the characteristic particle size relative to 
that responsible for the Galactic extinction law. This was 
done in previous papers (Goudfrooij et al. 1994 and Patil 
et al. 2007) by comparing the wavelength-dependence of the 
extinction efficiency of spherical dielectric grains with that 
of the observed R\ extinction values. This method relies on 
observations indicating that the Ra curves usually vary ap- 
proximately linearly with inverse wavelength, as observed 
for the Galactic extinction curve and for extinction curves 
of different galaxies as derived by Goudfrooij et al. 1994 and 
Patil et al. 2007. Therefore, the different extinction curves 
seem to run along or parallel to the Galactic curve. Consid- 
ering the optical part of the spectrum, and the upper and 
lower cutoff sizes estimated for the grains, Goudfrooij et al. 
(1994) assumed Q ex t on — (see also van de Hulst 1957). 



http 
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However, considering eq. |(HJ and eq. ((9j), it seems that this 
assumption leads to a situation where R\ is grain size in- 
dependent. Therefore, irrespective of the grain model used, 
all R\ values remain constant and higher order terms of the 
series expansion of the extinction efficiency are necessary in 
order to estimate the characteristic particle size relative to 
that responsible for the Galactic extinction curve. 

Adopting the MRN grain model, different extinction 
curves can be reproduced by evaluating the extinction ef- 
ficiency under the assumption of spherical grains and using 
the available Mie theory code derived from BHM1E (Bohren 
& Hoffman 19830. This fitting scheme should include sev- 
eral free parameters such as the size distribution power law, 
the abundance and the upper and lower size cutoff values 
(a+ and a-) for both graphite and silicate grains (see also 
Mathis & Wallenhorst 1981; Steenman and The 1989, 1991). 

Considering the extinction efficiency calculations and 
the different extinction curves, which we find to run paral- 
lel to the Galactic extinction curve, we can generally con- 
clude that increasing (decreasing) the number of upper end 
size grains with respect to the Galactic grain population 
produces larger (smaller) R\ values. Such an increase of 
the average grain size may be the consequence of a flat- 
ter size distribution power law, or of an increase in a+ and 
a_. Therefore, we conclude that extinction curves that lie 
above (below) the Galactic extinction curve correspond to 
larger (smaller) characteristic grain sizes with respect to the 
typical Galactic grain size. However, we note that varying 
the abundance ratio of the two silicate-graphite components 
may also account for different extinction values. Calculations 
A s 

show that increasing the — — ratio, as well as reducing the 
A g 

characteristic grain size, leads to a decrease of the extinction 
values and vice versa. 

We note that our measured extinction values in the 
near-IR to near-UV region cannot properly constraint all 
the free parameters mentioned above. Therefore, a precise 
fit should include a better determined extinction curve and 
an extension of this extinction curve to the UV region. Since 
the size distribution power law suggested by the MRN grain 
model appears to be consistent with the power law expected 
from grain-grain collisions (Biermann & Harwit 1980), and 
since the observed silica and carbon abundances in the ISM 
are approximatelly equal, we consider here only a variation 
of o+. As far as a_ is concerned, surface effects and quantum 
mechanics effects make Mie calculations suspect for particles 
smaller than about 0.01 /im. 

Using eq. |8j and the extinction cross-section calcula- 
tions we obtain the upper size cutoff value for the graphite- 
silicate grains mixture from the best fits to the extinction 
values for each sample galaxy and for the Galaxy. We find 
the average upper size cutoff value in our sample galaxies 
(see Table 2} to correspond only very approximately with 
the value of 0.22 /jm obtained for the MW. The a+ values 
are mostly lower than 0.22 pm, with only one (AM1118-290) 
significantly higher. 



4.3 Dust mass estimation 

In order to derive the dust mass we use the dust column 
density from 



f 4 3 

Sd = IdX / -ira p d n(a)da 



(11) 



where pd gives the specific grain density and Id represents 
the dust column length along the line of sight. The value of 
Id x riH can be inferred from the measured total extinction 
and the calculated extinction efficiency in the V band us- 
ing eq. ([9]l and a+ and a_ obtained for each sample galaxy. 
The specific grain densities of graphite and olivine-like 'as- 
tronomical silicate' grains are conservatively estimated as 
2.26 gr cm -3 and 3.3 gr cm -3 , respectively (DL84). The to- 
tal mass is derived by integrating the dust column density 
over the image areas occupied by dust lanes and is therefore 
given by Md — E^xArea, where M d is expressed in solar 
mass units. We also note that it is also possible to derive the 
optical dust mass using the alternative method described in 
§ 13.31 i.e., using the differences in the extinction and cross- 
sections in two optical bands. However, the latter requires 
subtracting the constant dust-free galaxy colours from the 
colour values where dust is present in order to obtain the 
correct difference of extinctions. 

The dust mass can be estimated independently from 
the far-infrared emission of each galaxy by using the dust 
grain temperature calculated from the IRAS flux densities at 

60pm and lOO^m using T d = (j^) (Young et al. 1989). 
The dust mass is then computed for each sample galaxy 
using the relation (Hildebrand 1983) 



(12) 



2 The code was written by Bruce T. Draine and is available at 
ftp: / /ftp. astro. princeton.edu/draine 



where a, p d and D are the grain radius, specific grain mass 
density and distance of the galaxy in Mpc (assuming Ho =70 
km s _1 Mpc -1 ), respectively; F u , Q„ and B v (Td) are the 
observed flux density, grain emissivity and the Planck func- 
tion for the temperature Td at frequency u, respectively. 
We can evaluate the graphite-silicate mixture mean emis- 
sivity by using published extinction efficiencies values and 
assuming thermal equilibrium. Since the quantity a/Q v is 
independent of a for a » A (Hilderbrand 1983), we calcu- 
late Q v for a mixture of 0.1pm graphite and silicate grains 
at A = lOO^im using the tabulated data from Draine (1985). 
These dust mass estimates represent lower limits, since IRAS 
was insensitive to dust cooler than about 20K which emits 
mostly at wavelengths longer than 100/xm. Table [J] lists the 
estimated dust mass from the total optical extinction and 
the estimated dust mass based on IRAS flux densities taken 
from the catalog of Knapp et al. (1989) for bright early- type 
galaxies. 



5 DISCUSSION 

We have studied the extragalactic dust extinction in nine 
dust- lane early- type galaxies. The Rv values derived above 
are close to the standard Galactic value, at least for most 
of the derived extinction curves. The mean value for these 
galaxies, Rv ~ 2.82 ±0.38, corresponds to the value derived 
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Table 4. Dust properties 



Object ow^C/mi) IRAS flux (mjy) Log Log(^f-) 

\ JVJ / d, optical \ / d. 







60/im 


lOO^m 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC3283 


0.19 ±0.01 


1200 ± 68 


4440 ± 207 


5.67 ±0.03 


6.48 ± 0.08 


NGC3497 


0.18 ± 0.01 


278 


< 1347 


6.03 ± 0.03 


< 6.43 


AM1118-290 


0.30 ±0.01 


190 ± 35 


800 ± 128 


6.26 ± 0.02 


6.76 ± 0.25 


NGC4370 


0.18 ± 0.01 


990 ± 40 


2900 ± 77 


4.46 ± 0.03 


4.96 ± 0.05 


AM1352-333 


0.16 ± 0.01 


110 ±35 


130 ± 365 


5.84 ±0.04 


4.54 ±2.83 


NGC5626 


0.20 ±0.01 


210 ±23 


760 ± 69 


5.93 ± 0.03 


6.44 ±0.15 


AM1459-722 


0.22 ± 0.01 






5.96 ± 0.03 




NGC5799 


0.10 ± 0.01 


320 ± 21 


2370 ± 294 


5.17 ±0.07 


6.96 ± 0.21 



by Patil et al. (2007) for galaxies with well-settled dust lanes, 
i.e., 2.80, supporting an explanation that the characteristic 
grain size responsible for the optical extinction in such galax- 
ies is slightly smaller than that in the MW, but within the 
error bars. These Rv values, with the exception of AM1118- 
290, are also in the range of values derived by Goudfooij et 
al. (1994), i.e., 2.1 to 3.3, for which they estimated a char- 
acteristic dust grain size up to 30% smaller compared to 
standard Galactic dust grains. We note that the extinction 
value derived for NGC5799 has an exceptionally low value 
of Rv = 1-90 ± 0.01; this is similar to the lowest extinction 
value obtained by Patil et al. (2007) and Goudfrooij et al. 
(1994). 

As described in § 13.21 the effects produced by some fore- 
ground light cannot account for lower Rv values and, if it 
all, dust uniformly embedded with stars or a dust lane in a 
distinct location within the galaxy will have the effect of pro- 
ducing larger values of Rv for the observed extinction law. 
Another potential pitfall in the derivation of the extinction 
values is the possibility of star formation in the dust lanes. 
Goudfrooij et al. (1994) discussed this effect thoroughly and 
concluded that, due to the light emitted from young massive 
stars, one may underestimate the extinction, especially at 
the shorter wavelengths. Thus, the real extinction curve may 
be somewhat more concave than the measured one and with 
larger extinction values. However, this may not account for 
the relatively low extinction values of NGC5799, for which 
Brosch et al. (1990) concluded from the relatively low ratio 
between the blue and 60/im flux densities that it does not 
form stars at a very high rate in its ISM torus, and therefore 
star formation should have a negligible effect on its derived 
extinction values. This may also be concluded for the rest of 
the sample galaxy. For most galaxies, we may also neglect 
the effect of forward scattering in dusty regions on the ob- 
served extinction curve. As already discussed in § 13.21 this 
effect leads to apparently lower extinction values towards the 
blue, making the extinction curve natter than the Galactic 
curve, whereas such flattening is not evident in any of our 
sample galaxies. 

We obtained extinction curves for eight of our nine sam- 
ple galaxies, with the exception of NGC5898. Concerning 
the latter, for which Patil et al. (2007) obtained Rv = 
3.15 ±0.46, the difficulty in determining its extinction curve 
may lie in the short SALT exposure times. We also found 

A 

that NGC3283 and AM1118-290 show noisy — - plots in 

Av 

most bands, suggesting a possibly complex dust structure 



or the presence of significant forward scattering. The lat- 
ter seems more reasonable for AM1118-290 which produces 
larger R\ values with respect to the Galactic extinction val- 
ues. For each of the NGC5799 and AM1459-722 galaxies 
Brosch et al. (1989) derived a flat extinction curve in the 
B & V filters, which did not correspond with the Galactic 
extinction trend or values. It is therefore reassuring to de- 
rive extinction curves which run parallel to the Galactic ex- 
tinction curve, though we find NGC5799 to have significally 
smaller grain sizes compared with the Galactic dust. The 
NGC5626 extinction law derived by Goudfooij et al. (1994), 
with Rv = 2.15, implies significally smaller grains in com- 
parison with Rv = 2.98 we found, the latter being closer to 
the Galactic value. The extinction values for NGC3479 are 
smaller than those derived by Patil et al. (2007) in the V 
and I bands, but still within the error bars. The major-axis 
dust lane elliptical NGC4370 shows an extinction curve very 
similar to that of the Galaxy. While the effect of forward 
scattering may be important in the case of major-axis dust 
lane disk galaxies, it does not seem significant for the case of 
pure ellipticals such as NGC4370. Such a possible effect, as 
well as that of a more complex dust structure, might reduce 
the observed total extinction of NGC4370 with respect to its 
optical depth. We did not find previously determined extinc- 
tion values for the NGC3283, AM1118-290, NGC4370 and 
AM1352-333 galaxies with which to compare our results. 

We introduced an alternative method to measure the 
extinction values by comparing colour-index maps. 

Although less informative, this method does not require fit- 
ting a model for the underlying galaxy, thus its results are 
not biased by the fitting procedures and may be compared 
with the traditional method. The comparison shows that the 

E(B V) vames are S enerau y consistent with the Galactic ex- 
tinction law and with the trends of the matching R\ curves. 
An exception is NGC3283, which shows a similar deviation 
from the Galactic extinction law in both extinction curves. 

Although we started this study with the intention of 
finding differences between the properties of extragalactic 
dust in our objects and the properties of MW dust, we are 
compelled by the results to conclude that no major differ- 
ences were found. It is possible that such differences, if they 
exist at all, may be detected in the space-UV. 

Using the optical extinction values we also evaluated 
the dust mass for each of the sample galaxy, and found it 
to lie in the range ~ 10 4 to lO 7 A/0. This is in good agree- 



Determining the extragalactic extinction law with SALT 13 



ment with the earlier estimates for early-type galaxies with 
dark lanes (Brosch et al. 1990; Goudfooij et al. 1994; Sahu 
et al. 1998; Ferrari et al. 1999; Dewangan et al. 1999; Tran 
et al. 2001; Patil et al. 2002 and Patil et al. 2007). The 
use of the optical method to derive dust masses is based on 
the assumption that the dust forms a foreground screen for 
the galaxy. This means that we do not consider the dust 
to intermix with stars within the galaxy, nor the possible 
twisting of dark lanes around or through the host galaxies. 
Therefore, this method provides only a lower limit to the 
true dust content of the host galaxies and it should be com- 
pared with studies using other methods. Comparison of dust 
masses derived using optical extinction with those derived 
using IRAS flux densities reveals that the latter are up to 
an order of magnitude higher than the first. An exception 
is AM1352-333 where the dust mass obtained using far-IR 
data is poorly determined. Our results are in good agreement 
with previous estimates (Goudfooij & de Jong 1995; Sahu 
et al. 1998; Dewangan et al. 1999; Patil et al. 2007), which 
also concluded that dust masses derived using IRAS flux 
densities are usually larger by up to an order of magnitude 
than those estimated using the optical method, so that the 
discrepancy is more significant for ellipticals than for lan- 
ticulars. Moreover, since IRAS is not sensitive to cold dust, 
which emits the bulk of its radiation longwards of 100 /tm, 
the IRAS dust mass estimates are a lower limit for the true 
dust mass which may even be an order of magnitude higher 
(Temi et al. 2004). Goudfrooij & de Jong (1995) suggested 
that this discprepancy may be solved if the interstellar dust 
exists in early-type galaxies as a two-component medium, 
where the less massive component is optically visible in the 
form of a dust lane or dust patches, and the more massive 
component is the diffuse dust which is distributed over the 
galaxy. 



6 CONCLUSIONS 

We presented SALT observations using the SALTICAM of 
early-type galaxies with dust lanes obtained in order to de- 
rive the extragalactic extinction law down to the ultravi- 
olet atmospheric cutoff. Our work also extends that done 
by previous authors since our derived extinction laws reach 
to A -1 w 2.94. This extension demonstrates the necessity of 
performing such studies at as short a wavelength as possible. 

Extinction curves derived for the sample galaxies run 
parallel to the canonical MW curve, implying similiar prop- 
erties between MW canonical grains and dust in the extra- 
galactic enviroment. We derived the ratio between the total 
V band extinction and the selective B and V extinction Rv 
for each galaxy, and obtained an avergage of 2.82 ± 0.38. 
These results are in agreement with Patil et al. (2007), in- 
dicating that galaxies with well-defind dust lanes tend to 
have slightly smaller Rv values with respect to the Galac- 
tic value, but within the error bars. This suggests that the 
characteristic grain sizes responsible for the optical extinc- 
tion are similar in size to the canonical MW grains. We veri- 
fied our results by using an alternative method to derive the 
extinction values, which is based on comparing colour-index 
maps. 

The dust content derived by the optical extinction 
method for the sample galaxies is in the range ~ 10 4 to 



10 7 Mq; this also is in agreement with Patil et al. (2007) 
and with previous studies. 
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Figure 6. Filled contour maps for NGC3283 with 0.5 magnitude steps. The images are 112" X 112" 
with north up and east to the left. The dark lane crosses the galactic center almost vertically. 




Figure 7. Filled contour maps for NGC3497 with 0.5 magnitude steps. The images are 112" X 112" 
where north is up and east is left. The dark lane is below and to the left of the galactic center. 
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Figure 8. Filled contour maps for AM1118-290 with 0.5 magnitude steps. The images are 112" xll2" 
where north is up and east is left. The dark lane is below and to the left of the galactic center. The 
surface magnitude level at 25 mag/arcscc 2 is labelled in the B-band image using values from the 
RC3 catalog. 




Figure 9. Filled contour maps for NGC4370 with 0.5 magnitude steps. The images are 112" X 112" 
where north is up and east is left. The dark lane crosses the galactic center almost horizontally. The 
surface magnitude level at 25 mag/arcsec 2 is labelled in the B-band image using values from the 
RC3 catalog. 
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Figure 10. Filled contour maps for AM1352-333 with 0.5 magnitude steps. The images are 
112" X 112" where north is up and east is left. The dark lane crosses above the galactic center hori- 
zontally. 




Figure 11. Filled contour maps for NGC5626 with 0.5 magnitude steps. The images are 112" xll2" 
where north is up and east is left. The dark lane is above and to the left of the galactic center. The 
surface magnitude level at 25 mag/arcsec 2 is labelled in the B-band image using values from the 
RC3 catalog. 
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Figure 12. Filled contour maps for AM1459-722 with 0.5 magnitude steps. The images are 
112" X 112" where north is up and east is left. The galaxy is centered in these images and the 
dark lane crosses horizontally below the galactic center from. The bright object above the galaxy is 
a nearby bright star. The surface magnitude level at 25 mag/arcsec 2 is labelled in the B-band image 
using values from the RC3 catalog. 




Figure 13. Filled contour maps for NGC5799 with 0.5 magnitude steps. The images are 112" xll2" 
where north is up and east is left. One dark lane is below the galactic center while another one is 
to the right of the galactic center. The surface magnitude level at 25 mag/arcsec 2 is labelled in the 
B-band image using values from the RC3 catalog. 
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Figure 14. Filled contour maps for NGC5898 with 0.5 magnitude steps. The images are 112" xll2" 
where north is up and east is left. The dark lane is to the right of the galactic center. 
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Figure 15. B-I colour-index maps of the sample galaxies. Brighter regions near the images center represent the 
suspected dust lanes. 



